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ABSTRACT 

Context. Only recently it has become possible to measure the thermal emission from hot-Jupiters at near-Infrared wavelengths using 
ground-based telescopes, by secondary eclipse observations. This allows the planet flux to be probed around the peak of its spectral 
energy distribution, which is vital for the understanding of its energy budget. 

Aims. The aim of the reported work is to measure the eclipse depth of the planet HAT-P-lb at 2.2/jm. This planet is an interesting case, 
since the amount of stellar irradiation it receives falls in between that of the two best studied systems (HD209458 and HD189733), 
and it has been suggested to have a weak thermal inversion layer. 

Methods. We have used the LIRIS instrument on the William Herschel Telescope (WHT) to observe the secondary eclipse of HAT- 
P-lb in the K,-band, as part of our Ground-based secondary eclipse (GROUSE) project. The observations were done in staring 
mode, while significantly defocusing the telescope to avoid saturation on the K=8.4 star. With an average cadence of 2.5 seconds, we 
collected 6520 frames during one night. 

Results. The eclipse is detected at the 4-cr level, the measured depth being 0.109+0.025%. The uncertainties are dominated by 
residual systematic effects, as estimated from different reduction/analysis procedures. The measured depth corresponds to a brightness 
temperature of 2136+}^K. This brightness temperature is significantly higher than those derived from longer wavelengths, making 
it difficult to fit all available data points with a plausible atmospheric model. However, it may be that we underestimate the true 
uncertainties of our measurements, since it is notoriously difficult to assign precise statistical significance to a result when systematic 
effects are important. 

Key words, techniques: photometric - stars: individual: HAT-P-1 - planetary systems 



1. Introduction 

Measurements of the secondary eclipse of an exoplanet, the mo- 
ment it passes behind its host star, allow us to probe the prop- 
erties of the atmosphere on the day-side of the planet. The first 
successful secondary eclipse measure ments have been obtained 
with the Spitzer Sp ace Telescope (Char bonneau et ail 2005: 
iDeming et alJ l2005). followed by many m ore secondary eclipse 
measurements from 3 .6^m to 2%m (e.g. iKnutson et al. 12008; 
iMachalek etai]l2008l) . see also the review bv lDemind(l2009h . 

These Spitzer observations indicate that so me planets ex- 
hibit a thermal inversion in their atmospheres (e.g. Knutso n~et al.l 
2008). This inversion is apparent when molecular bands in 
the infrared change from absorption to emission, resulting in 
a different shape of the planet's spectral ener gy distribution. 
iFortnev et al.l (120081) and iBurrows et al.l d2008) have proposed 
that such an inversion layer could be due to the presence of 
a strong optical absorber high in the planet's atmosphere, de- 
pending on the amount of stellar radiation the planet receives. 
Only for the planets receiving very high levels of irradiation, 
the hypothetical absorbing compound can stay in the gas-phase 
high up in the atmosphere, absorbing the stellar light very ef- 
ficientl y and causing the th ermal inversion. As a possible ab- 
sorber, Huben vet al.l J2003) suggested TiO and VO, although 



* Photometric time-series are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



recent work (Spiegel et al. 2009) show that it might be difficult 
to keep these molecules in the gas phase. At lower levels of stel- 
lar irradiation, the compound possibly condenses out and is sub- 
sequently removed from the higher layers of the atmosphere. 
Fortn ev et al.l d2008l) dub these two classes pM and pL respec- 
tively, in analogy with L and M stellar dwarfs. 

With the increase in the number of exoplanets studied with 
the Spitzer Space Telescope, it became clear that the above 
scheme cannot be solely dependent on the level of stellar ir- 
radiation. There are planets apparently exhibiting an inversion 
layer which receive less light from their host-star than required 
to keep the proposed absorber in the gas-phase (e.g. XO-lb, 
Machalek et al. 2008), while there are also planets that receive 
very high levels of irradiation, which d o not apparently ex - 
hibit an inversion layer (e.g. TrES-3b, iFressin et al.l 1201 61) . 
Recently, Kn utson et alj j2010l) . showed that there appears to be 
a trend between the presence of an inversion layer and the stellar 
activity (as determined from the calcium lines), with only plan- 
ets around stars with low acti vity exhibiting an inversion layer. 
Mad husudhan & Seagerl (1201 Ol) showed that the inference of a 
thermal inversion is not robust, and can also depend on the chem- 
ical composition of the atmosphere. 

Recently, several groups have presented measurements for a 
number of exoplanets in the near-infrared using ground-based 
telesc opes dde Mooii & Snellen! l2~009t ISing & Lopez-Morales! 
1 2009 1 ; iGillon et alJ 120091; iRogers et d] 120091; lAnderson et al l 
|20 lot lAlonso et al.ll2010l;lGibson et al.ll2010HCroll et all Eo 1 0a; 
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lLopez-Morales et all 120101: ICroU et al.1 l2010bi 120111) . Data in 
this wavelength region are very interesting because they probe 
the thermal emission at, or close to, the peak of the spectral en- 
ergy distribution of hot-Jupiters, which is vital for the under- 
standing of their energy budgets. 

Here we present the second result from the GROUnd-based 
Secondary Eclipse project (GROUSE), which aims to use tele- 
scopes at a variety of observatories for exoplanet secondary 
eclipse observations. As part of this survey we have already 
publi shed our detection of the secondary eclipse of TrES- 
3b (Ide Mooii & Snellenl2 009). which we retrospectively include 
as paper I. 

In this paper we present obser vations of the sec ondary 
eclipse of the hot- Jupiter HAT-P-lb dBakos et al.ll2007l) in K s - 
band as obtained with the Long -slit Intermediate Reso lution 
Infrared Spectrograph (LIRIS; lAcosta-Pulido et all 120021) in- 
strument on the William Herschel Telescope (WHT). HAT-P- 
lb, with a mass of Q.52Mj llp and a radius of l.2Rj up , orbits 
its GOV stellar host with a period of P=4.5 days at a dis- 
tance of 0.055 AU. This large orbital separation places HAT- 
P-l on the lower edg e of the proposed pM/pL transition bound- 
ary of Fortnev et al ] (120081) . Recent observations of this planet 
bv lTodorov et a l. (2010), who used the IRAC instrument on the 
Spitzer Space Telescope to determine the infrared brightness of 
HAT-P-1, showed evidence for a weak inversion layer. Its host 
star appears not to be very active dKnutson et alj|20loh . which 
provides an alternative explanation for the planet's inversion 
layer, in the case that stellar variability is one of the driving fac- 
tors for the structure of the atmosphere. 

In Sec. [2] we present our observations and data-reduction. In 
Sec. [3] we will present and discuss the results, and in Sec. |4] we 
will give our conclusions. 



2. Observations, data reduction and analysis 

The secondary eclipse of HAT-P-1 was observed on October 
2, 2009 using the imaging mode of the L ong-slit Intermediate 
Resol ution Infrared Spectrograph (LIRIS; lAcosta-Pulido et al] 
|2002|) instrument on the William Herschel Telescope (WHT) on 
La Palma. The expected time of mid-eclipse was 2l' ! :55'"UT. 
The observations started at 19'':59 m UT and lasted for 6.5 hours, 
although ~ 15 minutes were lost because the observer got locked 
up in the bathroom due to a broken door handle. Since this 
caused a gap in the photometry, we excluded the final 630 frames 
from further analysis. The exposures were taken in sequences of 
50 or 100 frames. The first two frames of every sequence are 
known to suffer from a reset-anomal)Q, visible by a very strong 
gradient in the background. We excluded the frames affected by 
the anomaly from further analysis. 

The LIRIS detector was windowed to one sub-array of 512 
by 512 pixels in order to increase the observing cadence. The 
pixel-scale of LIRIS is 0.25 arcsec per pixel, hence the field of 
view for these windowed observations was 128 by 128 arcsec - 
onds. Since HAT-P-lb has a slightly brighter companion star 
nearby, the field of view was enough to allow us to observe 
both HAT-P-1 (K=8.86) and the reference star, ADS 16402 A 
(K=8.41), simultaneously. 

The observations were performed in staring mode, keeping 
the stars as much as possible at the same position on the chip. 
Despite guiding, a slow drift of 3 pixels over 4 hours is visible, 
and subsequently a sudden jump of 5 pixels in right ascension 



1 http://www.ing.iac.es/astronomy/instruments/liris/liris_cookbook.pdf 
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Fig. 1. The flux in the highest pixel of the PSF of the brightest 
star as a function of time. 



occurs. The stars also slowly drifted by 2 pixels in declination 
during the night. 

The exposure times were varied between 1.3 and 2 seconds, 
in order to keep the flux-levels of the star in the linear regime 
of the detector. To allow these relatively long exposure times for 
these rather bright stars, we strongly defocused the telescope. 
This resulted in a donut-shaped PSF with a diameter of 15 pix- 
els (3.75 arcsec). Despite the defocus, the PSF was still small 
enough to avoid overlap between the target and the reference 
star (which are separated by ~ 1 1 .2 arcsec), although the diffrac- 
tion spikes originating from the support of the secondary mirror 
did overlap for certain periods (see section |2~4V 

2.1. Crosstalk, non-linearity corrections and flat-fielding 

As a first step in the data reduction we removed the intra- 
quadrant crosstalk, which is present as an additional flux in each 
pixel at the level of 10~ 5 of the total flux along its row. We subse- 
quently performed a non-linearity correction, for which we ob- 
tained a new set of measurements in August 2010. A plot of 
the value of the highest pixel in the PSF of the brightest star is 
shown in Fig. Q] For the non-linearity correction we observed 
several sequences of K r band dome-flats with varying exposure 
times to which we subsequently fitted a second order polynomial 
to the non-linearity for each pixel of the detector. Our current re- 
sults show that the detector is non-linear for all fluxes, reaching 
~8.5% at 30000 ADU, above this flux level the non-linearity in- 
creases rapidly. We used this non-linearity solution to correct all 
our frames. We excluded all frames with the peak flux of the 
brightest star above 30000 ADU. 

All the frames were flatfielded using a flatfield created from a 
set of dome flats, both with the dome-lights on and off to remove 
both the effects of emission from the telescope structure, as well 
as structure due to the dark-current from the detector. The fl at- 
field images were corrected for the intra-quadrant crosstalk and 
non-linearity in the same way as the science images. 

2.2. Removal of bad-pixels 

During most of the observations, one hot-pixel is located within 
the PSF of the target. This hot-pixel was corrected for by replac- 
ing it with the flux of the corresponding pixel in the reference 
star, scaled by the flux ratio of the stars. Since the position dif- 
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Fig. 2. The variation of the sky background during the night. 



ference between the two stars is not an integer, and the precise 
flux in inner edge of the donut is a strong function of position, 
we interpolated the PSF of the reference star to the same grid as 
HAT-P-1. before measuring the flux in the corresponding pixel. 
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Fig. 3. Sample of PSFs for our observations at high and low 
count-levels for different mount rotator angles. Left panels: 
mount rotator angle of -79.2 degrees. Right panels: +55.0 de- 
grees. The spokes caused by the support structure of the sec- 
ondary mirror can be clearly seen to rotate due to the alt-az 
mounting of the telescope. 



2.3. Background subtraction 

The background level during the observations varied, and the 
variation of the background per pixel is shown in Fig. A back- 
ground map was constructed from a separate set of observations 
of a different field taken on the same night. After subtracting this 
background map from the images, a small gradient was still ap- 
parent along the y-axis, which we believe is due to different ex- 
posure times (and dark currents) for the images used for the cre- 
ation of the background map and the science images. To remove 
this gradient, we first masked out all the stars and bad-pixels and 
subsequently determined the mean along a row of the detector, 
rejecting outliers deviating by more than 3 times the standard de- 
viation from the median of the row. The resultant profile along 
the y-axis was then fitted with a polynomial, in order to create 
a smooth background map, which was subsequently subtracted 
from all the columns of the image. 

The overall impact of the subtraction of a background image 
on the results is small, but not subtracting the background image 
results in an increased noise-level. 

2.4. Diffraction spokes from the secondary mirror support 

The close proximity of the reference star causes the diffraction 
pattern from the support structure for the secondary mirror to 
overlap with the other star's PSF. Since the WHT is on an alt-az 
mount, the instrument is rotated to maintain the same position 
angle on the sky. This causes the diffraction spikes at certain 
times to overlap, depending on the position-angle of the mount. 
In Fig. [3] two images are shown, on the left side with a rotator 
angle of -79.2 degrees, with the diffraction spikes not overlap- 
ing, while in the right side the rotator angle was +55 degrees, 
where the diffraction spikes do overlap. To account for this, the 
amount of flux in a region offset by the same amount as the sep- 
aration between the two stars was measured, and subsequently 
subtracted from the measured flux values of the stars, after ac- 
counting for the difference in flux levels. 
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Fig. 4. Data points show the uncorrected lightcurve of HAT-P- 
lb with the WHT on La Palma. The solid line indicates the best 
fitting model where we fit simultaneous for both the systematic 
effects as well as the eclipse depth. 



2.5. Aperture Photometry 

Subsequently, we performed aperture photometry on the two 
stars, with a radius of 14 pixels (3.5 arcsec), chosen to maximise 
the flux, but minimise the influence of the background. Any 
residual background was determined by measuring the mean flux 
in two boxes located 60 pixels (15 arcsec) above and below the 
midpoint between the two stars, with a box-size of is 131 by 41 
pixels (32.75 by 10.25 arcsec). In these boxes, we excluded bad- 
pixels from the background-determination, as well as masked 
the areas that can be affected by the diffraction spikes from the 
secondary mirror support. 
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Fig. 5. Different components fitted to the lightcurve due to the 
eclipse and systematic effects. Shown are, from top to bottom the 
eclipse model, the airmass, the exposure time and the x position 
on the detector. The different parameters are offset for clarity. 



Table 1. Parameters for HAT-P-1. References 
(l2007h: 



(2) iJohnson et alJ (|2008); 



(4) TodoroxetalJ^OlO) 



m lWinn etal.1 
(3) this work; 



parameter 


Value 


Reference 


Semi-major axis a (au) 


0.0551±0.0015 


(1) 


R, (R ) 


1.115+0.043 


(1) 


Rp (Rjup) 


1.204±0.051 


(1) 


i (deg) 


86.22±0.24 


(1) 


P (days) 


4.4652934±0.00072 


(2) 


b 


0.701±0.023 


(1) 


R„/R. 


0.11295+0.00073 


(2) 


Eclipse depth at 2.2/im 


0.109±0.025% 


(3) 


Eclipse depth at 3.6/im 


0.080+0.008% 


(4) 


Eclipse depth at 4.5/um 


0.135±0.022% 


(4) 


Eclipse depth at 5.8/um 


0.203+0.031% 


(4) 


Eclipse depth at 8.0/um 


0.238+0.040% 


(4) 
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Fig. 6. The lightcurve for the secondary eclipse of HAT-P-lb 
corrected for systematic effects. Overplotted is our best fitting 
eclipse model. 
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Fig. 7. The same as Fig. [6] but binned by 41 points. 



2.6. Correction for systematic effects 

The lightcurve of HAT-P-lb, corrected for the hot-pixel, diffrac- 
tion spikes, and normalised with that of the reference star, is 
shown in Fig. |4] Although the eclipse is already visible, the 
lightcurve is still clearly influenced by systematic effects, which 
we found to be correlated with the x-position on the detector, 
the exposure time and the geometric airmass. We fitted for these 
effects simultaneously with the eclipse depth. Other parameters, 
such as the background level, were tried, but resulted in no dif- 
ference in the final lightcurve. For the model of the se c ondar y 
eclipse, we used the formalism from Ma ndel & Agoll (|2002), 
where all the orbital parameters were kept fixed to the litera- 
ture values (given in Table [TJ, with the limbdarkening set to 0, 
and the timing of the secondary eclipse fix ed. We did not fit for 
a time offset, because iTodorov et al.l (12010) already showed that 
the timing of the secondary eclipses measured using the Spitzer 
Space Telescope shows no indication of any offset, using higher 
signal-to-noise data. 

We performed the fit both using the standard multi-linear re- 
gression algorithm in IDL (regress. pro) as well as using a Monte- 
Carlo Markov Chain (MCMC) analysis. After an initial fit, we 
performed a clipping of the residuals by excluding all points 
within bins of 38 pixels, for which the binned value deviate by 



more that 1.5-10" 3 . In this way we removed 228 unbinned points. 
After this clipping we performed an MCMC analysis. For the 
MCMC analysis, we created 5 chains with a length of 2 million 
steps each, using different starting positions. We trimmed away 
the first 200.000 points of each chain, to make sure that the initial 
starting point no longer influence the measured parameters. 

The best fitting model, overplotted on the data, is shown in 
Fig. [4] while the contributions for the individual components are 
shown in Fig. [5] 

3. Results 

The final, corrected lightcurve is shown in Fig. |6] with its 
binned (by 41 points) version in Fig. Q Assuming pure, white, 
Gaussian noise, we find an eclipse depth of 0. 1089+0.01 % (1 1- 
cr). However, clear residual systematic effects are visible in the 
binned lightcurve. The level of this red noise was estimated 
by comparing the standard deviation in the unbinned data with 
that of gradually increased binned lightcurves. From this we es- 
timate a red noise component at the level of 0.023%, which 
we added in quadrature to our eclipse depth uncertainty to be- 
come 0.109±0.025 %. We also pe rformed the residu al permua- 
tion ("Prayer Bead") method (e.g. lGillo n et al. 2007), where we 
add the best fit model to the residuals, after shifting the resid- 
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Fig. 8. The temperature-pressure (T-P) profile for the our at- 
mospheric models used for modelling the secondary eclipse 
depth. The black ( solid) line shows the uninverted model 
from iFortnev et alJ (I2008I) . the green (dashed-triple dotted) 
li ne is for a model wi th an ad-hoc temperature inversion as 
in iTodorov et al.l (l2010l) . the red (dashed) line shows the tem- 
perature pressure profile for an atmosphere with a normal tropo- 
sphere and an inversion, and the orange (dash-dotted) line is for 
an atmosphere with Venus-like clouds. The temperatures for the 
latter profiles are chosen to provide a reasonable fit to the IRAC 
bands longward of 4//m 



uals by Appoints. Points are wrapped around when performing 
the shift. By refitting the data for all possible shifts between 
and the number of data points, and measuring the best-fit eclipse 
depth, we can get an estimate for the true uncertainties in the 
data. We find a depth of 0.1089+°°^ %, which is fully consis- 
tent with the uncertainties estimated from the red-noise. 

This results overall in an eclipse depth of 0.109+0.025 %, 
which corresponds to a brightness temperature in the K s -band 
of T/,=2136^qK. This brightness temperature is significantly 
higher than the expected equilibrium temperature of this planet, 
which lies between T eq = 1540 K for an albedo of and inefficient 
energy redistribution from the dayside to the nightside (P„=0), 
and 1 eq -\\9Q K for an albedo of 0.3 and efficient redistribu- 
tion of energy from the dayside to the nightside (P„=0.5). The 
higher brightness temperature could be due to the fact that in the 
Kj-band we are looking deeper into the atmosphere, where the 
layers are warmer. 

ITodorov et al.l d2010l) observed the secondary eclipse of 
HAT-P-lb with the IRAC instrument on board the Spitzer 
Space Telescope, at 3.6, 4.5, 5.8 and 8.0 /urn, and find eclipse 
depths of 0.080+0.008%, 0.135+0.022%, 0.203+0.031% and 
0.238+0.040% in the different bands respectively. Combining 
our own measurement of the secondary eclipse in K s -band with 
the IRAC measurements, we can construct a spectral energy 
distribution (SED) from the near-infrared to the mid-infrared, 
which is shown in Fig. [9] and compare it to atmospheric models. 



3.1. Atmospheric models 

We modelled the thermal emission spectra for a few ad hoc cases 
to try to understand what type of atmosphere would be needed 
to fit the measured eclipse depths, including our new measure- 
ment in Ks-band. We used a radiative transfer model that can cal- 
culate multiple scattering of the thermal radiation, since we in- 
cluded a case with scattering clouds (see below). The scattering 
is calculated using the doubling-adding method, which makes 
use of the fact that one can calculate the reflection and trans- 
mission properties of a combination of two atmospheric layers 
from the properties of the individual layers. For each layer in the 
model atmosphere we started with an optically thin layer with 
the same scattering properties as the model layer. For the opti- 
cal thin case analytical expressions are available to calculate the 
layer transmission and reflection. The thin layer was then dou- 
bled several times to match the real optical depth of the model 
layer. Subsequently, the different model layers were added to 
calculate the spectra for the entire i nhomogeneous atmos phere. 
For more explanations we refer to IWauben et al.l ([1994), who 
describe the model in great detail. Spectra were calculated at 10 
different emission angles, from which a disc-averaged spectrum 
was calculated. 

We included absorption of water and CO, whose ab- 
sorption properties were tak en from the (old) HITEMP 
database (Rothman et al.lll995h . Although the water data in this 
database is not the most up-to-date, a line-b y-line comparison in 
K-band with the new HITEMP database iRothman et alj [2010) 
yielded only differences many times smaller than the errors on 
the measured eclipse depths, which we neglect for this qualita- 
tive assessment. For speedy calculations, we pre-computed ab- 
sorption properties at a ra nge of temperatures an d pressure using 
the correlated-k method (lLacis & Oinasl Il99lh with a spectral 
resolution of 0.01 fim between 1.8 /vm and 5 /mi and 0.1 fim be- 
tween 5 yum and 10/vm, and a Voigt line shape. Assumed mixing 
ratios of water and CO were 2e-4 and 5e-4 respectively for all 
cases. 

The models extend from 1.8jum to lOjum, and were made 
for different temperature-pressure (T-P) profiles and composi- 
tions. _WeJiav^jrmdejnodels for the T-P profile for this planet 
from Fortnev et al.l (120081) . which does not include an inversion 
layer, as well as for models with an inversion layer. We also var- 
ied the composition of the atmosphere and made a model with 
Venus-like clouds. In Fig. [8] we show the temperature-pressure 
profiles for the different models, and in Fig. [9] we show the ex- 
pected eclipse depths within the observed bands for the different 
models, binned to a resolution of 0. 1/im. 

Fig. |9] shows that an a tmosphere calculated from first prin- 
ciples dFortnev et al.l 120081) cannot reproduce th e Spitzer mea- 
surem ents, and that a weak inversion is needed (Todorovet al. 
2010). However, a weak inversion throughout the atmosphere 
produces a very low planet signal in Ks-band, unlike what we 
derive from our measurements. The reason for this is that Ks- 
band is a region of very low absorption and hence we see deep 
in the planet atmosphere, where temperatures are relatively low 
in the inversion-only case (green lines). Contribution functions 
for the different wavelengths are shown in Fig. [10] To increase 
the eclipse depth in Ks-band the lower atmosphere needs to be 
made hotter. However, this also increases the eclipse depth in 
L-band, another spectral region of low absorption (red cureve 
in Fig. [9j. L-band emission can be suppressed somewhat by in- 
cluding large amounts of CQ 2 and CH4 (calculate d here from 
HITEMP and HITRAN 2008 dRothman et all 20091) data respec- 
tively), but this also reduces the signal in Ks-band (blue curve in 
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Fig. 9. Models for the emission from the atmosphere of HAT-P-lb, for the T-P profiles of Fig. [8] using the same colour scheme. 
The diamonds indicate the expected flux within the observed bands for the different models. The points with errorbars indicate the 
observational data. The dotted lines indicate the transmission curves of the different bands. In addition, a blue (long dashed) curve 
is shown for an atmosphere with the same T-P profile as the red (short dashed) curve, but with a high concentration of methane and 
carbon-dioxide ( 1 ■ 10 4 and 2-10 -5 respectively). 
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Fig. 10. Normalised contribution functions for the atmospheric emission of HAT-P-lb in the K s -band and the different IRAC chan- 
nels, using the same colour scheme as in Figs. [8] and [9] although the different linestyles are now used for different bands. The 
distri bution functions are normalised by their sum. The left panel shows the contribution functions for the non-inverted model 
from lFortnev et al.l (120081) . the middle panel shows the contribution functions for the model s with both a t roposp here and an inver- 
sion layer, and the right panel shows the same for the ad-hoc temperature inversion used in lTodorov et al.l(l2010l) 
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Fig. [9]). Reducing the water or CO abundances also do not pro- 
duce a better fit. Hence, we cannot fit all data points well using 
any clear atmosphere with gases that are predicted to be abun- 
dant on hot exoplanets. One remedy to fit both the Spitzer data 
and our Ks-band eclipse depth is to include a Venus-like cloud. 
The clouds on Venus are made of concentrated sulphuric acid 
droplets, which have the property that they are strongly scatter- 
ing below wavelengths of ~3/mi and strongly absorbing above 3 
^m, with only little va riation of extinction with wavelength (e.g. 
Grinsnoo n et al.llT993h . Because of the scattering nature of the 
clouds, radiation from Venus' hot lower atmosphere still reaches 
space at the night side below 3fim, despite an optically thick 
cloud layer surrounding the planet. Above 3 /mi the night side 
emission originates from the clouds. To show this potential effect 
for HAT-P-lb, we inserted a cloud layer with an optical depth of 
1 .5 at the tropopause with properties identical to those of Venus' 
1-yum sized cloud-particles. Indeed, brightness temperatures at 
Ks-band are higher than anywhere else in Fig. [9] and all data 
points could potentially be fitted well if the lower atmosphere 
is made even hotter. However, we could not find a physically 
plausible candidate for cloud materials that could mimic Venus' 
clouds on a hot exoplanet. So, at present we do not find any 
suitable atmosphere that could explain the high Ks-band eclipse 
depth. 
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4. Conclusion 

Using the LIRIS infrared camera on the WHT, we determine the 
eclipse depth of the extrasolar planet HAT-P-lb in K s -band to be 
0.109+0.025% (~4<r), with the uncertainties in the eclipse depth 
dominated by residual systematic effects, as estimated from dif- 
ferent reduction/analysis procedures. The measured depth cor- 
responds to a brightness temperature of 2136^[^jK. This bright- 
ness temperature is significantly higher than those derived from 
longer wavelengths, making it difficult to fit all available data 
points with a plausible atmospheric model. It may be that we 
underestimate the true uncertainties of our measurements, since 
it is notoriously difficult to assign precise statistical significance 
to a result when systematic effects are important. 
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